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r\ • Abstract 

en 

Gluon fusion is the main production mechanism for Higgs bosons with masses 
CN ' up to several hundred GeV in pp collisions at the CERN Large Hadron Col- 

J^ ' lider. We investigate the effects of the CP-vioIating phases on the fusion 

lO ■ process including both the sfermion-loop contributions and the one-loop in- 

duced CP-violating scalar-pseudoscalar mixing in the minimal supersym- 
metric standard model. With a universal trilinear parameter assumed, every 
^^ I physical observable involves only the sum of the phases of the universal trilin- 

■^ ' ear parameter A and the higgsino mass parameter //. The phase affects the 

Q.' lightest Higgs-boson production rate significantly through the neutral Higgs- 

JL , boson mixing and, for the masses around the lightest stop-pair threshold, 

D i it also changes the production rate of the heavy Higgs bosons significantly 



through both the stop and sbottom loops and the neutral Higgs-boson mix- 
ing. 
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I. INTRODUCTION 

The experimental observation of Higgs bosons and the detailed confirmation of their 
fundamental properties are crucial for our understanding of the mechanism responsible for 
the electroweak symmetry breaking and they constitute one of the most important exper- 
imental programme at the CERN Large Hadron Collider (LHC), the physics potential of 
which has been investigated with a considerable amount of effort [0]. Within the minimal 
super symmetric standard model (MSSM), the lightest Higgs boson is expected to be well 
within the reach of the LHC and also the heavier Higgs bosons with their masses less than 
1 TeV could be found at the LHC. Gluon fusion [Q is the main production mechanism for 
Higgs bosons with masses up to several hundred GeV at the LHC. The dominant produc- 
tion of these three neutral Higgs bosons, which proceeds through loops of both quarks and 
squarks (primarily, those of top and bottom flavor), can be significantly affected by non-zero 
CP-violating phases in the MSSM [0,||, associated with the higgsino mass parameter fj, and 
trilinear scalar couplings Af. 

It has recently been pointed out ^ that these CP-violating phases do not have to 
be suppressed in order to satisfy the constraints from electron and neutron electric dipole 
moments (EDMs). In this context, the so-called effective supersymmetry model is of 
particular interest, where sfermions of the first and second generations are decoupled, but 
sfermions of the third generation remain relatively light. Based on the scenarios of this type, 
many important works in B decays M , dark-matter searches, and collider experiments 



have been recently reported. In the present work, we concentrate on the effects of the CP- 
violating phases on the Higgs-boson production in gluon-gluon fusion at the LHC including 



both the sfermion-loop contributions |^ and the neutral Higgs-boson mixing |Tl . 

In this scheme, the characteristic CP-violating phenomena in the Higgs sector are a 
possibly large mixing between the CP-even and CP-odd neutral Higgs bosons and an 
induced relative phase ^ P,p!I| between the vacuum expectation values of the two Higgs 
doublets. In particular, the heavy CP-even and CP-odd Higgs bosons experience a maximal 
mixing in most cases due to their sharp mass degeneracy, in particular for large masses. In 
this light, the neutral Higgs-boson mixing cannot be simply neglected with the argument 
that it is induced at the one-loop level. Neglecting the other small loop contributions 
from the chargino and neutralino sectors, the Higgs potential of the MSSM has one-loop 
radiative corrections mainly from the stop and sbottom sectors due to the large t and b 
Yukawa couplings. 

On the other hand, the contributions of the sfermion loops to the gluon-gluon-Higgs 
vertices are expected to be large for sfermion masses similar to or smaller than the Higgs- 
boson mass. Concerning these contributions, Dedes and Moretti ^ have shown that the 
squark loop content of the dominant Higgs production mechanism via the gluon-gluon fusion 
mode at the LHC could be responsible for large corrections to the known cross sections. 
However, the possible large modifications due to the neutral Higgs-boson mixing are not 
discussed in their work. 

We re-visit the important production mechanism of Higgs bosons via gluon-gluon fusion 
at hadron colliders including both the stop/sbottom loop diagrams and the CP-violating 
mixing among three neutral Higgs bosons, and provide a detailed analysis of the effects 
of the CP-violating phases on the Higgs-boson production rates. The one-loop Feynman 



diagrams contributing to the process gg -^ Hi in the MSSM are shown in Fig. 1. 

From the viewpoint of radiative corrections, the neutral Higgs-boson mixing contributes 
to the Higgs-boson production in gluon-gluon annihilation at the two-loop level. For the 
sake of consistency, it might be necessary to take into account all the two-loop vertex 
corrections to the gluon-gluon-Higgs interactions. However, the most dominant radiative 
corrections are expected to be largely contained in the propagators and the vertex corrections 
are much smaller in size compared to the propagator corrections |T^. This point has been 
demonstrated explicitly in the Higgs-boson pair production process e~^e~ — > HiHj by Demir 
| |TB| . Moreover, the sharp degeneracy of the heavy CP-odd and CP-even Higgs states causes 
a large mixing between two heavy Higgs-boson states even for a tiny CP-violating correc- 
tion to the Higgs-boson mass matrix. In this light, we could neglect the two-loop vertex 
corrections to the gluon-gluon-Higgs interactions to a good approximation and concentrate 
on the effects of the neutral Higgs-boson mixing modifying the Higgs-boson propagators. 
Clearly, it will be important to demonstrate quantitatively through explicit calculations that 
the two-loop vertex corrections are negligible compared to the neutral Higgs-boson mixing 
in the Higgs-boson propagators, even if we do not touch on the comparison in the present 
work. 

The paper is organized as follows. In Sec. II we describe our theoretical frameworks; the 
MSSM Lagrangian, the neutral Higgs-boson mixing and the sfermion mixing, including a 
brief discussion of experimental limits on the parameters of the model. Sect. Ill is devoted to 
the general description of the Higgs production via gluon-gluon fusion at hadron colliders. 
In Sec. IV we present our numerical results. And then we summarize our findings and draw 
conclusions in Sec. V. Finally, the Appendix is devoted to listing all the Feynman rules 
relevant for our numerical analysis. 



II. CP VIOLATION IN THE MSSM 



A. The MSSM Lagrangian 



In this section, we list the relevant MSSM Lagrangian terms including CP-violating 
sources, and the effective Lagrangian for the Higgs potential with radiative interactions to 
quartic couplings coming from enhanced Yukawa couplings of the third generation. For a 
general purpose, it is useful to replace two conventional MSSM Higgs doublet fields Hi and 
H2 by $1 = +iT2Hl and $2 following the notations used in Ref. |11]. Then, the relevant 
MSSM Lagrangian is divided into four parts: 

- Aoft = rhlQ^Q + mllJ*lJ + mlb*b + (/i,A<I>1QD - htAt^lir^QU + H.c.) , (1) 

-Cf = \h^W + \ht^2^^Q\^ - {i^KQ^^2b* + fxh;QHT2^u* + h.c.) 

-(/i^D*$fir2 + h;U*<l>l){hbiT2<l>lD - ht<l>2U) , (2) 






2|$f ^TsQp + 2|$^Q|2 - gtg($t$^ + $t$^^ 



,/2 



+ '—m^2 



$I$i; 



liQ^Q)-^iU*U) + ^{D*D) 



-Cf = hh 



&ii(tL,&L)'$I+H.C. +ht tR(tL,6L)'ir2<l>2 + H.C. 



(3) 
(4) 



with the notations Q^ = [tLjbi)'^, U* = tn, and D* = h^ for the scalar top and bottom 
quarks. 

The MSSM introduces several new parameters in the theory that are absent from the 
SM and could, in principle, possess many CP-odd phases |^. Specifically, the new CP 
phases may come from the following parameters: (i) the higgsino mass parameter /z, which 
involves the bilinear mixing of the two Higgs chiral superfields in the superpotential; (ii) 
the soft SUSY-breaking gaugino masses M^ (a = 1,2,3), where the index a stands for the 
gauge groups U(l)y, SU(2)i and SU(3)c, respectively; (iii) the soft bilinear Higgs mixing 
masses 77^12 5 which is sometimes denoted as -B/i in the literature; and (iv) the soft trilinear 
Yukawa couplings Aj of the Higgs particles to scalar fermions. If the universality condition 
is imposed on all gaugino masses at the unification scale Mx, the M^'s have a common 
phase. Likewise, the different trilinear couplings Af are all equal at Mx-, i.e. Aj = A. 

The conformal-invariant part of the supersymmetric Lagrangian has two global U(l) 
symmetries; the U(1)q Peccei-Quinn symmetry and the V{1)r symmetry acting on the 
Grassmann-valued coordinates. As a consequence, not all phases of the four complex pa- 
rameters {/i, m^25 ^a, A} turn out to be physical, i.e. two phases may be removed by redefin- 
ing the fields accordingly. Employing the two global symmetries, one of the Higgs doublets 
and the gaugino fields can be re-phased such that Ma and 771^2 become real. In this case, 
arg(/i) and arg(y4) are the only physical CP-violating phases at low energies in the MSSM 
supplemented by universal boundary conditions. 

We work in the theoretical framework provided by the MSSM, including explicitly the 
CP-violating phases $^ = arg(/i) and $^ = a.Tg{A) of the universal trilinear parameter. 
Moreover we take the case where the first and second sfermions are so heavy as to avoid 
the stringent EDM constraints and assume the universal soft-breaking mass parameters and 
phases for the third generation squarks: 

\At\ = \Ab\=A, $^^ = $^^ = <|.^. (5) 

Consequently, all the physical quantities are determined by {tan/5, |y4|, \fi\, $^, $^}, and two 
sfermion masses as well as the charged Higgs-boson mass mH±, as shown later. 



B. Neutral Higgs boson mixing 

The most general CP-violating Higgs potential of the MSSM can conveniently be de- 



scribed by the effective Lagrangian |11 



+Ai($l$i)2 + Aal'I'^a'^z)' + A3($I$i)($^$2) + A4($I$2)($|*i) 

-1-A5($1$2)' + Xli^l'^lf + A6($I$l)($I$2) + A*(<l'I<l>i)($^$i) 

+Xr{4^2){^\^2) + A;($^$2)($^2<^i) . (6) 

At the tree level, the kinematic mass parameters are given by /if = —mf — |/ip and /i2 = 
—ml — |yup, and the quartic couplings determined by the SM gauge couplings as 



Ai = A2 = -i(/+/), A3 = -i(/-9«) 



A. = -^s' 



A.c; — Afi — A7 



0, 



(7) 



^2y 



and m?9 are soft-SUSY-breaking parameters related to the Higgs sector. 



n2 



Here, mj, m 

Beyond the Born approximation, the quartic couphngs A5, Ag, A7 receive significant radiative 
corrections from trihnear Yukawa couphngs of the Higgs fields to scalar-top and scalar- 
bottom quarks. These parameters are in general complex and the analytic expressions of 
the quartic couplings can be found in the Appendix of Ref 



n 



In the existence of the CP-violating quartic couplings, the three neutral, physical Higgs 
bosons may strongly mix with one another. To describe this CP-violating mixing, it is 
necessary to determine the ground state of the Higgs potential to determine physical Higgs 
states and their self-interactions. To this end we introduce the linear decompositions of the 
Higgs fields 



$1 




+ iai, 



$9 



J( 



71' 



V2 + (f)2+ 10.2) 



(8) 



with t>i and f 2 the moduli of the vacuum expectation values (VEVs) of the Higgs doublets 
and their relative phase ^. These VEVs and the relative phase can be determined by the 
minimization conditions on Cy- It is always guaranteed that one combination of the CP-odd 
Higgs fields oi and 02 (G° = cos/5ai + sin/3a2) defines a fiat direction in the Higgs potential 
and it is absorbed as the longitudinal component of the Z field. As a result, there exist one 
charged Higgs state and three neutral Higgs states that are mixed in the presence of CP 
violation in the Higgs sector. Denoting the remaining CP-odd state a = — sin/?ai + cos/3a2, 
the neutral Higgs-boson mass matrix describing the mixing between CP-even and CP-odd 
fields is given by 



M" 



N 



Ml Ml 



MIp 



PS 



M% 



(9) 



where M^ and A^| describe the CP-preserving transitions a —* a and (</ 
respectively, and M^g = (M'^p)'^ contains the CP-violating mixing a 
analytic form of the sub-matrices is given by 



^i,<A'2; 



h,92) 



The 



Ml 



M'sp 



Mi 



m„ 



n{mi2e'^)+v 



SpCp 



2n{\^e''^)spcp + -7^(A6e^^)4 + -n{\je'^)s 



i^\ 



M\ 



m„ 



-SpCii 



-SpCp 
C/3 



2Ai4 + 27^(A5e2^«)s^ + 27^(A6e^«)s^C;3 A34S/3C^ + n{\e^y% + 7^(A7e^«)s^ 
A34S/3C^ + n{\S)cl + ll[\^S)sl 2A24 + 27^(A5e2^«)4 + 27^(A7e 



Correspondingly, the charged Higgs-boson mass is given by 



i€^ 



S/3C/3 



• (10) 



m^± =ml + ]^Kv^ - niX.e^'^y . (11) 

Taking this very last relation and the universality of the trilinear parameters into account, we 
can express the neutral Higgs-boson masses as functions of {mH±, |/i|, 1^41, $^, $^}, the com- 
mon SUSY scale MgusY, the parameter tan/3, and the physical induced phase ^. However, 
with the chargino and neutralino contributions neglected, the radiatively induced phase ^ 
accompanies the the /i parameter without exception so that only the sum of the three phases, 
^Afi = ^A + ^fi + ^, enters the loop-corrected Higgs potential. 

The CP-violating induced relative phase ^ between two Higgs doublets can be obtained 
analytically by combining the two relations for the minimization 0. First of all, we make 
use of the fact that a U(1)pq rotation allows us to take rri^2 to be real and for a notational 
convenience define Xq, 6, and 6'; 

Xq = Aq c^ + Ay s^ , 
6 = 




6'=^-^-^-X,jsin2f3. (12) 

Then, the CP-violating induced phase ^ is determined by the relations; 

sine = -^ [n{6)lCXe) - I{6)^\6\'-I'{Xe 

co«^ = +^ {^(5)^(A6) + nS)^/\6\'-I'{Xe)] , (13) 

and the soft-breaking positive bilinear mass squared 771^2 is given by 



9 V 

m 



12 
2 



|j(M')J(A6) + 7^(M') Vl^P - X2(A6) - \SmK)] . (14) 



The induced phase C, involves only arg(y4yu), i.e. $ = $^ + $^ and so vanishes when 
$ vanishes, even if each of the CP-violating phases might not have to vanish. On the 
other hand, the size of 6 or 6' is proportional to the pseudoscalar mass m^ to a very good 
approximation so that if large, i.e, decoupled, the induced phase ^ is diminished. Since 
the size of the induced phase is also inversely proportional to sin 2/3, the phase grows with 
increasing tan/3. 

Clearly, the CP-even and CP-odd states mix unless all the imaginary parts of the 
parameters A5, Ag, A7 are vanishing. Since the Higgs — boson mass matrix A4% describing 
the mixing is symmetric, we can diagonalize it by means of an orthogonal rotation O; 

0^M\0 = diag(m^3 , m^^ , m^J , (15) 

with the ordering mH^ < 1TLH2 ^ itT'Hs ■ The size of the resulting CP-violating neutral Higgs- 
boson mixing is characterized by the factor (l/327r^)(|/i||Aj|y^/MgugY) sin$^j.^, where Yf 
is the Yukawa coupling of the fermion /, ^Af^l = Arg(y4j/i) + ^ for / = t, b, and Msusy is a 
typical SUSY-breaking scale, which might be taken to be the average of the two sfermion 
(/i,/2) masses squared. 

6 



C. Sfermion mixing 

The sfermion mass matrix squared of the left- /right-handed sfermions is given by 

where, for the case of scalar top quarks, 

X 
At = m}^ +m1 + -(4m^ - m|) cos 2/3 , 

Bt = mt\A* -/ie*^cot/3| , 

2 
Ct = m}^ + m] + -m\sl^ cos 2(3 , 

0i = arg(A-/i*e"*^cot/5), (17) 

and, for the case of scalar bottom quarks, 

Ab = ml +ml- - (2m^ + m|) cos 2j3 , 

Bb = mh\Al — /ie*^ tan/9| , 

Cb = m-^^ + m^ - ^rrizS^, cos 2/3 , 

0b = arg(Ab-/i*e~*^tan/3), (18) 

where mf and rrij are the left /right-handed soft-SUSY-breaking stop/sbottom quark 

masses squared, respectively. The Hermitian mass matrix squared 7W| can be diagonalized 
by an unitary transformation 

Uj M)Uf = diag(mj-^ , m| ) [ mj^ < "^/J , ( 19) 

with the parameterization for the unitary matrix 

^ / cos^; -sin^;e-^/\ 
^ l^sin^/e*"^/ cos% ) ' ^ ^ 

satisfying — 7r/2 < ^/ < 0, and the sfermion mass eigenvalues and mixing angles are given 
by 

^/i.2 = ^susY T 2 A/ , 
2Bf 



sin 2^ 



/ 



{Af - CfY + AB} 



cos 2^. = ^ ^ , (21) 

l{Af -CfY + 4B} 



with the phenomenological parameters Mg^gy and Aj defined as 



^^-'susY — 2 ~ 2 




A/ = ml - ml = ^{Aj - Cff + AB} . (22) 

We note that there exists a freedom of redefining the sfermion fields by phase rotations. 
So, in constructing the Feynman rules for the Higgs-sfermion-sfermion vertices, we make 
the transformation of the sfermion basis (/l, fn) into the mass basis (/i, /2) by use of the 
unitary matrix Uf and a 2 x 2 phase rotation matrix Pf-. 

= P/f//(|), (23) 

where the 2x2 matrix Pf is defined in terms of the phase of the trilinear term Af as 

P^ = diag(e-'*■^//^ e'*-^//') . (24) 

The additional phase rotation renders the Higgs-sfermion-sfermion couplings dependent 
only on the sum of the phases, ^Arfi = ^Af + *^^ + ^- Furthermore, the orthogonal matrix O 
transforming the Higgs boson weak basis into the Higgs boson mass basis involves only the 
sum of the phases as noted before. As a result, with the universal trilinear term assumed, 
only one combination ^a/i = ^A + ^fi + ^ of the CP- violating phases dictates all the physical 
quantities related to the sfermion sectors and the CP-violating induced neutral Higgs-boson 
mixing. 

Not only the phase ^a^ but also the magnitudes of |A| and |/i| affect the values of the 
sfermion masses by modifying the mass splitting. Since we are interested mainly in the 
effects of the CP-violating phases, we take the real parameters to be fixed. In this case, 
the allowed range for the CP-phase $a/x is determined by the values of the real parameters 
through the constraint: 

A < MlusY • (25) 

In principle, A can be as large as 2MgugY) but the expressions of loop corrections to the 
Higgs potential |TI[] is valid when the tighter constraint ( |2^ ) is satisfied. With a simple 



assumption that two diagonal elements of the stop mass matrix. At and Ct, are equal, the 
constraint (|^) leads to the inequality for the phase ^Af^'- 



\A\' + \f^\y t^n' P - % 



SUSY 



showing that the CP-odd phase $yi^ could be constrained by the measurements of the scalar 
top-quark masses, depending on the values of real parameters \A\ and |/i|. 

Fig. 2 exhibits the allowed maximal value of the CP-violating phase $^^ on the plane 
of {\A\, |/i|/tan/?} for the scalar top-quark sector for Msusy = 500 GeV and m^ = 175 
GeV. A few comments are in order: (i) the full range of the phase is not allowed unless \A\ 
or |//|/tan/3 is less than M|ugY/2mt ^ 710 GeV; (ii) the upper-left and lower-right regions 
are not allowed due to the self-evident constraint | cos$^^| < 1; (iii) the lower-left hatched 



triangular area allows the full range for the phase $a/x; (iv) if At is not equal to Ct, the phase 
is more constrained. In addition, the fact that every physical quantity is determined by only 
the combined phase $a/x = ^a + ^^i + ^ implies that large CP-violating effects are caused 
when \A\ and |/i|/tan/3 are similar as well as large in size. One natural consequence is then 



that for a large tan/3, which is disfavored by the two-loop EDM limits []n[, the effects of 
the CP phases are reduced. Based on these observations, we consider two sets of {\A\, \fi\} 
for a fixed value of tan (3 = 3 and Msusy = 500 GeV in the following numerical analysis Q : 

RRl : {\fi\ =0.5TeV,|v4| = 0.5TeV}, 

RR2: {|/i| = 1.2TeV,|A| = 0.4TeV}. (27) 

We note that the full range for the phase ^a/i is allowed in the RRl, but the phase is 
restricted to the intervals [0°, 140°] and [220°, 360°] in the RR2. The Higgs-boson production 
rates, which we are interested in, are CP-even observables and so they are dependent on 
only the cosine of the phase ^Afj.- Keeping in mind this point we have considered only the 
range [0°, 180°] for the phase $^^ in Fig. 2. 

Experimentally, the tightest limits on the squark masses come from direct searches at 
the Tevatron. Concerning the ti mass, the limit on m^^ can be drawn at around 120 GeV or 
so for a large value of tan f3 assuming the mass of the lightest neutralino is smaller than 50 



GeV ||T6|. A less model-dependent mass bound of 90 GeV is obtained at LEP2 [|T^. As for 



the lightest sbottom mass, m^ , this is excluded for somewhat lower values [Q. Since both 
scenarios give the lightest stop mass larger than 200 GeV, the present experimental bounds 
on the sfermion masses can be safely neglected in our analysis. 

III. HIGGS BOSON PRODUCTION AT HADRON COLLIDERS 

A. Parton level production cross section 

The basic production processes of supersymmetric Higgs particles at hadron colliders are 
essentially the same as in the Standard Model (SM). Important differences are nevertheless 
generated by the modified couplings, the extended particle spectrum, and the existence of 
Higgs bosons of mixed parity. There are three crucial contributions to the Higgs-boson 
production in gluon-gluon collisions; (i) heavy fermion contributions, (ii) heavy sfermion 
contributions affected by sfermion mixing due to the large Yukawa couplings, and (iii) large 
Higgs-boson mixing due to the same large Yukawa couplings. The contributions (i) and 
(ii) are at the one-loop order, while the last contribution (iii) is at the two-loop order in 
character. However, the one-loop Higgs-boson mixing cannot be neglected as a higher-loop 
effect, because even a small loop contribution to the scalar-pseudoscalar mixing term can 



^ For large values of At — ///tan/3 and moderate m^ , T{gg — > H) can be siginificantly reduced 
by even an order of magnitude due to the destructive stop contributions |l^. This destructive 
interference effect is not so large in our choice of SUSY parameters. 



cause a maximal mixing between the heavier CP-even and CP-odd neutral Higgs fields due 
to their intrinsic close degeneracy in the MSSM. 

The leading-order (LO) partonic cross section aLo{99 ~^ Hi) [i = 1,2,3] for the gluon 
fusion of Higgs particles in the MSSM with explicit CP violation can be expressed in the 
narrow-width approximation by their scalar and pseudoscalar couplings, gl^j- and gpjj- to 
each fermion /, and their couplings g^ff, to each sfermion /: 



//' 



(^Loigg 



C^LO 




m 



H, 



4 - mj 



^fy'if 



+ 



E 



9pff 
ruf 



^pA^f) 



(28) 
(29) 



where the explicit form of the couplings {(7*^^,5'^^^, (?*;;} is given in the Appendix and the 
form factors A^/, Apj and A 7 can be expressed in terms of the scaling function f{Tix = 



Asfir) 



A. 



pjy 



r-'f{r), 



Af{r) 



r-'[-l+r-'f{r)], 
where the scaling function /(r) stands for the integrated function 



fir) 



I -^ln[l-Ary{l-y)] 



2 Jo y 



ln(l 



-VT- 

-vr- 



— ITC 



r<l, 
r > 1. 



(30) 



(31) 



For small tan/3 the contribution of the top loop is dominant, while for large tan/3 the bottom 
loop is strongly enhanced. The contributions of the squark loops can be significant for the 
Higgs-boson masses above the squark-pair thresholds [|19|. 

The limits of both large and small loop masses are interesting for SUSY Higgs particles. 
The contribution of the top loop to the Higg coupling can be calculated approximately in 
the limit of large loop masses, while the bottom contributions to the Higg couplings can 
be calculated in the approximation of small b masses. The limits of large loop masses are 
approximated with the following limiting values of the form factors. 



A. 



f 



2/3, 



A. 



pf 



A 



f 



1/3, 



(32) 



and the Apf is not altered by QCD radiative corrections due to the non-renormalization 
of the axial anomaly. Incidentally, glj-jr and gpj-j- are always proportional to the Yukawa 
coupling Yf, i.e. the fermion mass rrif so that the fermionic contributions remain finite 
even for large fermion masses. However, the sfermion-loop contributions vanish because of 
the suppression by the factor 1/m^. On the other hand, in the opposite limit where the 
quark-loop mass is much smaller than the Higgs masses, the amplitudes are the same for 
the scalar and pseudoscalar form factors: 



A 



sf 



A 



pf 



1 

'4^ 



ln(4r) + in 



(33) 



10 



This result follows from the restoration of chiral symmetry in the limit of vanishing quark 
masses. 

If the mass of the produced Higgs boson is smaller than the ti and tit^ pair thresholds, the 
form factors are real; above the thresholds they are complex. Fig. 3 shows the dependence 
of the real (solid) and imaginary (dashed) parts of the form factors, Agf (left), Apf (middle), 
and A? (right) on the Drell-Yan variable r. All the real parts are peaked at r = 1 above 
which all the imaginary parts become finite. Since the fermionic form factors are determined 
by the top quark mass rrit = 175 GeV, the production cross sections are expected to be 
enhanced near mn = 'ifnt with unsuppressed couplings of the Higgs boson to top quarks. 
Similarly, the sfermionic form factor will be enhanced for uih = 2mi with unsuppressed 
couplings of the Higgs boson to scalar top quarks. 

The Higgs bosons and sparticles of the MSSM that enter the process gg — > Hi can be 
also produced via other channels, both as real and virtual objects. From their search at past 
and present colliders, several limits on their masses and couplings have been drawn. As for 
the neutral Higgs bosons of the MSSM, the most stringent bounds come from LEP in the 
CP-invariant theories. For both mho and tti^o these are set at around 80 GeV by all the 



LEP Collaborations |J20| for tan/5 > 1. There is also a more recent and somewhat higher 



limit on m^jo from ALEPH |^ of about 85 GeV Q for tan/5 > 1 at 95% confidence level, 
using data collected at the cm. energies in the range between 192 GeV and 196 GeV and 
a total luminosity of about 100 pb^^. However, the existence of CP-violating phases could 
weaken the present experimental bounds on the lightest Higgs mass up to about 60 GeV 



11 



B. Hadronic level production cross section 

In the CP-preserving limit, the cross section near the top-quark-pair threshold are not 
reliable taking into account the perturbative two-loop QCD corrections which are infinite 



at the threshold |jT9|. 



It is well known |[19|J23|| that next-to-leading order (NLO) corrections to the processes 
gg — i> Hi from ordinary QCD are very large. However, it has been shown that they affect 
the quark and squark contributions very similarly [l^, and can be approximated by the so- 



called K factors, defined by the ratios of the higher-order cross sections to the leading order 
cross sections. They vary little with the masses of the neutral Higgs bosons in general, yet 
they depend strongly on tan/3. For small tan/?, their size is between 1.5 and 1.7; for large 
tan/? they are in general close to unity, except when the lightest Higgs boson approaches 
the SM domain. Therefore, we may take K to he a. constant for a given tan/3 to a good 
approximation 0. 



The most recent limit obtained by the LEP collaborations at the end of the '99 run which is 



still preliminary and not combined is 90 GeV |22|. 



Since the QCD correction to the CP-odd part of the production amplitude is singular at the 



top~quark-pair threshold, the cross sections are not realiable near the threshold |24|. 
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In order to find the physical cross section it is necessary to integrate the parton-level 
production cross section with the distribution of gluons in a proton, 



a 



{pp -^ Hi) = K dxidx2 g{xi, Q)g{x2, Q) o-lo{99 -^ Hi 



Apgg 
= KcrLoi99^H,)r^, (34) 

where g{x, Q) is the distribution of gluons in the proton and the the gluon-gluon luminosity 
is defined as 

T-— = J dx-g{x,Q)g{T/x,Q), (35) 

and the Drell-Yan variable r = mjj./s with s being the invariant hadron collider energy 
squared. We note that the cross section decreases with increasing Higgs mass. This is, 
to a large extent, a consequence of the sharply-falling gg luminosity for large Higgs-boson 
masses. Therefore, the significant modifications of the lightest Higgs mass by the radiative 
corrections cannot be neglected. 

The dependence of the cross sections on the Higgs-boson masses is encoded in the gluon- 
gluon luminosity function and the parton-level cross sections a{gg —>■ Hi), the latter of which 
is expected to have a weaker dependence. On the other hand, the lightest Higgs-boson mass 
is less than 130 GeV irrespective of the CP-violating phases and its lower limit can be set by 
its present collider searches. This limit is of course dependent on the CP-violating phases. 
Nevertheless, in our analysis we simply take 70 GeV for the lower limit of the lightest Higgs- 
boson mass. Fig. 4 shows the gluon-gluon luminosity function r(d£^^/dr) with respect to 
the mass of the produced Higgs mass at the Tevatron withy^ = 2 TeV (solid line) and at 
the LHC with a/s = 14 TeV (dashed line), taking CTEQ4m p5[ as the parton distribution 
routine. Note that the luminosity function is larger at the LHC at least by two orders of 
magnitude than at the Tevatron, and the difference is enlarged for large Higgs masses. So, 
the Higgs-boson production rates are much larger at the LHC and remain considerable for 
the Higgs-boson masses up to several hundred GeV, while the production rates are strongly 
suppressed at the Tevatron for large Higgs masses. 



IV. NUMERICAL RESULTS 

In our numerical analysis, we include the contributions from only the loops of the t, b, 
ii^2, and 61^25 that are indeed the dominant terms because of their large Yukawa couplings. 
We take tan/3 to be 3 and assume for the scalar top/bottom quarks masses MgusY = 500 
GeV and A = (500 GeV)^ for which the masses of the lighter and heavier scalar top quarks 
are approximately given by 

rrii^ = 350 GeV , m^-^ =610 GeV , (36) 

respectively. Therefore, the Higgs-boson production is dominated by the contributions from 
the top quark and the lighter scalar top-quark loops if the Higgs-boson masses are less than 
ITeV. 
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As shown in Sec. II, the CP-violating phases in the sfermion sectors affect the production 
rates of MSSM Higgs-bosons in two ways; through the direct loop contributions and through 
the neutral Higgs-boson mixing. In order to estimate the relative size of two contributions, 
firstly we consider the ratios of the production cross sections with the top-quark and sfermion 
loop corrections to those with only the top-quark loop contributions without including the 
neutral Higgs-boson mixing O: 

where the superscript indicates that the neutral Higgs-boson mixing is not included. In 
Fig. 5 we present the ratios with respect to the neutral Higgs-boson masses. The upper 
three frames are in the scenario RRl, while the lower three frames are in the scenario RR2. 
The phase $ = $a + $^ is taken to be 0° (sohd), 60° (dashed), 120° (dotted), and 180° 
(dot-dashed). We note in passing that in RR2 the allowed maximal angle of the phase is 
about 140°. Concerning the ratios (|37|) we make a few comments: 

(a) The stop-loop diagrams contribute little to the production rates unless the Higgs- 
boson masses are similar to or larger than the stop-pair threshold. 

(b) In the CP-invariant theories, H2 is the CP-odd pseudoscalar, which does not couple 
to the diagonal pairs of sfermions. So, the ratio is a unity for $ = irrespective of the 
Higgs-boson mass rri/^j as shown with the solid lines in the middle frames. 

(c) For H2, the production rates are enhanced due to the appearance of the H2tii\ coupling 
for the non-trivial CP-violating phase $ and the size of the enhancement is the 
same for $ = 60" and 120°, reflecting the fact that the diagonal H2iit\ coupling is 
proportional to sin 2$. 

(d) The CP-violating phase $ changes the ratio for H^ more signiflcantly than that of H2 
above the stop-pair threshold. This is because the new stop-loop diagram interferes 
with the top-loop diagram at the amplitude level only for the CP-even if 3, but not 
for the CP-odd H2. 

(e) A large value of |/i| enhances the sensitivity of the production rates of the heavy 
Higgs-bosons to the phase $ for rriH^^ > 2mi^, in particular, for the heaviest Higgs 
boson. 

Consequently, the stop-loop contributions can manifest themselves in the Higgs-boson pro- 
duction in gluon-gluon fusion only if the Higgs-boson mass is similar to or larger than the 
stop-pair threshold. 

Secondly, the neutral Higgs-boson mixing signiflcantly modifles the couplings of the 
Higgs bosons to fermions as well as gauge bosons in the presence of the CP-violating phases 



111 . Therefore, it is interesting to consider the ratios of the Higgs-boson production cross 
sections with the neutral Higgs-boson mixing to those without the mixing, while the top 
loop contributions but not the sfermion loop corrections are included: 
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Fig. 6 shows the ratios (|38D with respect to the produced Higgs-boson masses. The upper 
three frames are in the scenario RRl, while the lower three frames are in the scenario RR2. 
The phase $ = $a + $m is taken to be 0° (sohd), 60° (dashed), 120° (dotted), and 180° 
(dot-dashed), and in RR2 the phase angle is allowed up to about 140°. 
The ratios (^) in Fig. 6 present several interesting aspects: 

(a) In both RRl and RR2 the production rates for the lightest Higgs boson changes 
significantly with the phase $, reflecting the fact that the couplings of the Higgs 
boson to the top-quark pair is modified significantly through the neutral Higgs — boson 
mixing. 

(b) On the contrary, the production rates for the heavy Higgs bosons do not change so 
much except for the Higgs-boson masses around the top-quark pair threshold. It 
is also worthwhile to note that if the ratio for H2 is enhanced, the ratio for H2 is 
suppressed, and vice versa. This is a reflection of the fact that for a large Higgs-boson 
mass, the heavy Higgs bosons shows a two-state mixing. 

Consequently, if only the top-quark loops are considered, the neutral Higgs-boson mixing 
affects the production of the lightest Higgs boson significantly, but it does not affect the 
production of the heavy Higgs boson significantly unless the masses are around the top- 
quark pair threshold. 

Thirdly, in order to estimate the effects of the neutral Higgs-boson mixing to the Higgs- 
boson production with both the fermion and sfermion loops, let us consider the ratios of 
the full production cross sections with the neutral Higgs-boson mixing to those without the 
neutral Higgs-boson mixing, which have been considered by Dedes and Moretti [Q: 

^idlM , (39) 

^Uf + f] 

In Fig. 7 we present the ratios (^) of the Higgs-boson production cross sections with the 
neutral Higgs-boson mixing to those without the mixing with respect to the mass of the 
produced Higgs boson. The upper three frames are in the scenario RRl, while the lower 
three frames are in the scenario RR2. The phase $ = $a + ^^ is taken to be 0° (solid), 60° 
(dashed), 120° (dotted), and 180° (dot-dashed). We note that in RR2 the phase angle is 
allowed only up to 140°. 

Comparing Fig. 7 with Figs. 5 and 6 we observe several interesting aspects: 

(a) On the whole, the ratios are very sensitive to the CP-violating phase $ for large 1^41 
and |/i|, in particular, for a large value of |/i|. 

(b) The production rate of the lightest Higgs boson is modified mainly by the neutral 
Higgs-boson mixing as clearly seen by comparing Fig. 6 and Fig. 7. 
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(c) In RRl, the production rates for the heavier Higgs bosons are not affected by the 
neutral Higgs-boson mixing. However, they change a lot in the scenario RR2 with a 
large value of \fi\. The production rates are most sensitive to the CP-violating phase 
$ if Mh2 3 > 2mj^ ^ 710 GeV, i.e. the threshold for a light stop-pair. 

(d) Comparing Figs. 5 and 6 with Fig. 7, one can easily find that the couplings of the heavy 
Higgs bosons to sfermions also are significantly modified by the neutral Higgs-boson 
mixing, resulting in a large enhancement of (TLoidQ -^ H2) and a simultaneous large 
suppression of o-Lo{.gg -^ H^) for non-trivial values of the phase $. 

To summarize, although the neutral Higgs-boson mixing affects the Higgs-boson production 
by gluon-gluon fusion at the two-loop level, their contributions to the production rates can 
be very significant for some non-trivial values of the phase $, when the higgsino mass 
parameter \fi\ and the trilinear term \A\ are sizable. In addition, we emphasize that the 
neutral Higgs-boson mixing modifies the mass spectrum of neutral Higgs bosons as well as 
the couplings of the Higgs bosons to fermions, gauge bosons, and scalar fermions significantly. 

Finally, we consider the parton-level Higgs-boson production cross sections aLoi.99 ~^ 
Hi) (i = 1, 2, 3) including the effects of the CP-violating phases in both the sfermion loops 
and the neutral Higgs-boson mixing. In Fig. 8 we present the cross sections with respect to 
the Higgs-boson masses, mu^ (left), mu^ (middle), and mH^ (right) in RRl (upper) and RR2 
(lower) for four values of $; 0*^ (solid), 60" (dashed), 120° (dotted), and 180° (dot-dashed). 

Concerning the dependence of CFLo{,9g ^' H^) on the mass m^j and the phase $, we find 
several interesting aspects: 

(a) cTLoigg ~^ Hi) increases sharply with increasing mui because of the sharp increase 
of Asf with niHi for mu^ < 2mt ^ 350 GeV. In RR2, the CP-violating phase could 
suppress the cross section by a factor of five around itlhi = 90 GeV in spite of the 
enhancement by the neutral Higgs-boson mixing so that the possibility of finding the 
lightest Higgs-boson at the LHC might be seriously reduced. 

(b) cFLo{gg -^ H2), where H2 is a CP-odd pseudoscalar in the CP-invariant limit, is 
maximal at rriHi = 2mt, because Apf has a sharp peak at r = 1, i.e. the tt threshold. 

(c) CTLoigg -^ H2) is strongly modified by the CP-violating phase $ near the tit*i thresh- 
old. This is because a non-trivial coupling of the CP-odd content of H2 to diagonal 
stop pairs is developed and H2 is no longer a pure CP eigenstate due to the neutral 
Higgs-boson mixing in the existence of the phase $. There is also a little suppression 
in the production rates at the top-pair threshold due to the suppression of the coupling 
of iJ2 to top quarks. 

(d) cTLoigg ~^ H^) has a relatively dull peak at the ti threshold, as expected from the 
dependence of Asf on the Drell-Yan variable r in Fig. 4. 

(e) The CP-violating phase $ changes (yLo{gg ~^ H^) significantly near the tit* threshold; 
the enhancement is due to the stop-loop contributions, even though the effects are 
reduced by the neutral Higgs-boson mixing. 
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To recapitulate, the CP-violating phase $ changes (JLoidQ ~^ Hi) significantly over the 
whole range of the mass mij^ between 70 GeV and 110 GeV. d^Loiog ~^ H2) is maximal 
at the ti threshold, and its size is enhanced significantly near the threshold of the light 
stop-pair by the phase $. And, o-lo^QQ ~^ H^) is modified significantly at both the ti and 
titl thresholds. 

CONCLUSIONS 

We have investigated the effects of the CP-violating phases on the gluon-fusion process 
for the production of three neutral Higgs bosons including both the sfermion-loop (as well 
as fermion-loop) contributions and the one-loop induced CP-violating scalar-pseudoscalar 
mixing in the minimal supersjTumetric standard model with explicit CP violation. 

We have observed two generic aspects for CP violation from the stop/sbottom sectors: 

(a) If a universal trilinear parameter \A\ is assumed, the CP-violating phases $a and ^^ 
enter every physical quantity with only the sum $ = ^a + '^^• 

(b) The (measured) stop masses could constrain the CP-violating phase, depending on 
the values of \A\ and |/i|/tan/5. 

From a detailed numerical analysis, we have found that the CP-violating phase $ could 
reduce the production rate of Hi significantly, (which is dominated by the top-quark loop) 
and the phase enhances the production rate of H2 and H^ significantly near and above the 
stop-pair thresholds. Clearly, if the light stop mass is smaller, the effects of the CP-violating 
phase $ could be observed for smaller Higgs-boson masses. 

The most crucial observation of the present work is that the CP-violating neutral Higgs- 
boson mixing cannot be neglected, but should be included together with the third-generation 
sfermion loop contributions for reliable estimates of the Higgs-boson production rates in the 
presence of the non-trivial CP-violating phases. 

We conclude by noting that similar phenomena appear in the production of the Higgs 
bosons in two-photon fusion |26,27| as well as the two-photon decays of the Higgs boson. 



In this case, high energy colliding beams of linearly polarized photon beams, which can be 
generated by Compton back-scattering of linearly polarized laser light on electron/positron 



bunches of e~^e linear colliders |28], are expected to provide a clear-cut, direct means for 



determining the the CP properties of the produced Higgs bosons. We will discuss these 



crucial issues in detail in a future publication |29 . 
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APPENDIX A: FEYNMAN RULES 

In this Appendix, we list all the Feynman rules needed for the present work on the 
Higgs-boson production via gluon-gluon fusion. 

Firstly, the interactions of a gluon to fermions are given by 

Cgff = -gsfiY{T^)ijfjG;^, (A.l) 

and those of gluons to sfermions by 

C^jj = -igs{T'^)ijf*j d, fuG''^ + g',{T''T')jjf*jf,jG';,G^' , (A.2) 

where T" = X°'/2 (a = lto8) are the eight SU(3) generators in the adjoint representaion, 
and G'^ the corresponding eight gluon fields. 



Secondly, the interactions of the neutral Higgs-boson fields with fermions are described 
by the Lagrangian 

^^// = - n """"^o/ f K - ^^^475] fH^^ Mff + il,glff)f H, , (A.3) 

Obviously, the Higgs-fermion-fermion couplings are significant for the third-generation 
fermions, t, h and r. We readily see that the CP-violating neutral Higgs-boson mixing 
induces a simultaneous coupling oi Hi {i = \, 2, 3) to CP-even and CP-odd fermion bilin- 
ears // and fi'j^f- This can lead to a sizeable phenomena of CP violation in the Higgs 
decays into polarized top-quark or tau-lepton pairs. 

Thirdly, the Feynman rules for the Higgs-sfermion-sfermion vertices, involving all the 
mixing and phases and including the phase rotations of the sfermion fields, can be written 
in terms of C^.p^^ as 

^Hj,f\ = 9}.fJjfkHi, 

where a denotes three neutral Higgs-boson fields {a, 0i, ^2} and {j3, 7} denotes the chiralities 
{L,R}. The chiral couplings Ci-a^ can be obtained in a rather tedious but straightforward 
way and, for the scalar-top quarks, the couplings are given by 

CUl = , 

P^t ( r^^ \^ 

^a;RL ~ \^a;LR) ' 

CUr = , (A.6) 

^<j>i;RL = V^<f>r,LR) i 

i 

'3 



~4 _ gmw 2 2 .. „s 



gml ^ gmw ^ fl 2 2 



^<f>2;LL- ^ ^ "^ ^2 "^/3 I 9 o-^W 
mwSjS C^y \l 6 

~t _ gmt I ^ I 
"^f-^^^- 2mwsr'\' 

^<I>2;RL — \y'<f>2;LR) ■. 
,2 



gm^ gmw 2 2 
mwsp cfy 3 
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and, for the scalar-bottom quarks, they are given by 



^a;RL — \^a;LR) i 

Cl,RR = , (A.9) 

^b _ 9^1 gmw / 1 , 1 



^b _ 9jnb.. I 

P~ib I (^^ \* 

^4>r,RL — \^(f>Y\LR) ) 

~, _ gmj ^^^ 1 2 . , .(.^ 

'-'<)>r,RR- ^ ^ + ^2 C/3 S^y, (^A.iUJ 



gmw /I 1 2 
2;ii ~ ^2 -^/S I 2 3"^^ 



^4>2;RL — \^4,2;LR) ; 



i>2;RR ~ ^2 '^/3q' 






with $At^ = "^At + "^A* + 'C s-iid $Ab^ = *^Ab + $^ + ^. Therefore, if the univesal phase 
(^At = ^At is taken, the induced phase C, depends only on the combination ^a + $^ and as a 
result all the Higgs-sfermion-sfermion couplings are determined by only one CP-odd phase 
$ = $yi + $^, the sum of the two phases $^ and <l>^. 
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FIG. 1. The one-loop diagrams contributing to the process gg ^ Hi] (A) SM-hke contributions 
from fermions /, in particular, top and bottom quarks, and (B) SUSY-hke contributions from 
sfermions /, in particular, scalar top and bottom quarks in the MSSM. 
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FIG. 2. The allowed regime for the CP-violating phase $a^ = ^a + ^^i + S, on the parameter 
space of two real SUSY parameters \A\ and |;u|/tan/3 for a fixed value of tan/3 = 3. The lower-left 
hatched triangular area allows the full range for the phase $a^) while the lower-right and upper-left 
regime are not allowed due to the obvious condition | cos^a^iI ^ 1- The number of each contour 
line denotes the maximally allowed value of the phase for given \A\ and |^|/tan/3. 
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FIG. 3. The real (solid) and imaginary (dashed) parts of the form factors, Asf (left), Apf 
(middle), and A? (right), with respect to the Drell-Yan variable r. The vertical dotted line is 
to denote the thresholds for Higgs-boson decays into a fermion-pair or a sfermion-pair. The real 
parts have their peaks at r = 1 and the imaginary parts are non-zero only above the threshold. 
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FIG. 4. The gluon-gluon luminosity function with respect to the Higgs-boson mass niH at the 
Tevatron (sohd) with ^/s = 2 TeV and at the LHC (dashed) with ^ = 14 TeV. The Drell-Yan 
variable r denotes mjj/s by using the CTEQ4m parameterization for gluon densities. The lumi- 
nosity function is larger at the LHC by two orders of magnitude than at the Tevtaron. It decreases 
very abruptly with increasing Higgs-boson mass niH at the Tevatron, while it does not decrease 
so much at the LHC. 
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RRl : l|il=0.5 TeV , IAI=0.5 TeV 
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RR2 : l|il=1.2 TeV , IAI=0.4 TeV 
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FIG. 5. o'^oif + f]/^Lo{f] with respect to the mass of the produced Higgs boson. The upper 
three frames are in the scenario RRl, while the lower three frames are in the scenario RR2. The 
phase $ = $A + <^M is taken to be 0° (soUd), 60° (dashed), 120° (dotted), and 180° (dot-dashed). 
Note that in RR2 the phase angle only up to about 140° is not allowed. 
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RRl : l|il=0.5 TeV , IAI=0.5 TeV 
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FIG. 6. <5"Lo [/] /'5'io [/] '^ith respect to the mass of the produced Higgs boson. The upper three 
frames are in the scenario RRl, while the lower three frames are in the scenario RR2. The phase 
$ = $^ + $^ is taken to be 0° (solid), 60° (dashed), 120° (dotted), and 180° (dot-dashed). 
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RRl : l|il=0.5 TeV , IAI=0.5 TeV 
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RR2 : l|il=1.2 TeV , IAI=0.4 TeV 
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FIG. 7. aLo[f + f]/^Lo[f + /] with respect to the mass of the produced Higgs boson. The 
upper three frames are in the scenarios RRl, while the lower three frames are in the scenario 
RR2. The phase $ = ^^ + $/, is taken to be 0° (sohd), 60° (dashed), 120° (dotted), and 180° 
(dot-dashed) . 
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FIG. 8. CTLo{gg ~^ Hi) {i = 1,2,3) with the neutral Higgs-boson mixing as well as the 
stop/sbottom loop corrections with respect to the mass of the produced Higgs boson. The upper 
three frames are in the scenarios RRl, while the lower three frames are in the scenario RR2. The 
phase $ = $A + <I'^ is taken to be 0° (sohd), 60° (dashed), 120° (dotted), and 180° (dot-dashed). 
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